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GZK horizon and magnetic fields 

Etienne Parizot 
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15, rue Georges Clemenceau, 91406 Orsay Cedex, France 

We discuss some aspects of the propagation of high-energy cosmic rays (CRs) in turbulent magnetic fields, and 
propose a formula for the diffusion coefficient based on accurate simulations in a wide energy range. We discuss the 
transition between ballistic and diffusive regimes and the magnetic horizon of CR sources, showing how magnetic 
fields of a few nG could modify the shape of the GZK feature. Such fields would roughly be in equipartition with 
the extragalactic CRs, and could build through the resonant growth of waves within the lifetime of the universe. 



1. Introduction 

Two independent problems must be distin- 
guished in the physics and phenomenology of 
ultra-high-energy cosmic rays (UHECRs). The 
first one relates to their production and involves 
the identification of efficient astrophysical or ex- 
otic sources. The main issues are the prediction 
of a (sufhciently high) maximum energy, a defi- 
nite source spectrum and source composition, and 
global energetics. Most models do not currently 
have such a wide predictive power. The main 
challenge for astrophysical models appears to be 
the acceleration of particles up to energies around 
3 10^° eV, while the so-called top-down models 
encounter some difficulty to "hide" other types 
of radiation associated with the production and 
propagation of UHECRs. 

The second, independent problem relates to 
UHECRs propagation in the universe, which in- 
volves two complementary aspects: energy losses 
(mostly via interactions with the CMB) and par- 
ticle deflection in ambient magnetic fields. In the 
case of nuclei, a third issue is related to compo- 
sition changes, as photo-nuclear interactions lead 
to the "erosion" of primary heavy nuclei. 

As described forty years ago by Greisen^ and 
by Zatsepin and Kuzmin|5] (GZK), a sudden de- 
crease in the observed flux of cosmic-rays is ex- 
pected when the energy of the CMB photons in 
the rest frame of the propagating proton (or nu- 
cleus) reaches values in excess of the pion pho- 
toproduction (or photodisintegration) threshold. 



This should occur around 8 10^^ eV. Whether 
such a "cutoff" has been observed or not is still 
unclear, as the two main experiments which could 
address this question, AGASA and HiRes, do not 
agree on the answer, although with a statistical 
significance of only ~ 2(7 [21 . The determination 
of the CR spectrum above ^ 10^" eV is one of the 
major goals of the Pierre Auger Observatory 0], 
which should provide its first results soon. 

2. UHECRs and magnetic fields 

In the original GZK argument and in most 
subsequent discussions, the extragalactic mag- 
netic field has been assumed to be very low, say 
B <, 10~^^ G, so that its effect on particle tra- 
jectories could be neglected. Such an assumption 
is likely to be correct, indeed, as no mechanism 
has been proposed to produce higher magnetic 
fields over very large scales. Current data sug- 
gest that, except inside galaxies and possibly in 
the internal parts of galaxy clusters - which rep- 
resent a negligible fraction of the volume of the 
universe - magnetic fields are generally well below 
the equipartition values which would be required 
to produce significant deflections. However, it is 
expected that the sources of UHECRs are located 
in high magnetic field regions, because virtually 
all energetic phenomena take place in strongly 
magnetized regions, and intense magnetic fields 
are certainly required for particle acceleration. 
Therefore, even if they are globally negligible in 
the universe, magnetic fields may play a role in 
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the phenomenology of UHECRs, not only at ac- 
celeration, but also during their transport, and it 
is interesting to ask how they modify the origi- 
nal GZK argument and in particular change the 
shape of the GZK feature at ultra high energy. 

An important remark should also be made here. 
It is well known that CRs accelerated in super- 
nova remnants by the diffusive shock acceleration 
mechanism produce themselves the MHD waves 
thanks to which they diffuse back and forth and 
cross the shock front many times. The growth of 
Alfven waves by the so-called streaming instabil- 
ity occurs whenever an anisotropic flux of CRs is 
found with a global streaming velocity larger than 
the Alfven speed. Just as the CRs are scattered 
by resonant waves having a wavelength compa- 
rable with their gyroradius, the waves which are 
amplified by a CR stream have wavelengths of the 
order of their gyroradius in the underlying field. 
In a simplified approach, one may evaluate the 
growth rate of the waves (in the linear regime) 
by writing that the decrease of the CR momen- 
tum parallel to the stream must be compensated 
by the increase of the wave momentum. 

During a scattering time Ts ~ the reso- 

nant CRs are deflected by ~ 90°, which amounts 
to a momentum transfer of scrv jc? per unit vol- 
ume, where £cR is the CR energy density. The 
wave energy and momentum densities are related 
by = PwVA, where wa = Bo/^/Awp is the 
Alfven velocity in the regular field, Bq. Thus: 
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If Bl is the wave field, so that — B\/2p,Q, 
and Tg = p/qBq is the gyroradius of the resonant 
particles, we have 
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and d-Uw/di ~ Uw(wA/rg)(ecR/uo) (for v 
The wave growth timescale is thus: 
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where we have introduced the "equilibrium" mag- 
netic field, -Bcq, which may be reached if the waves 



can grow up to a point when the magnetic en- 
ergy density is similar to that of the cosmic rays 
(at a given resonant energy). Modeling the CR 
distribution function (in the energy range of in- 
terest) as $(£') ~ 3 lO^^eV^m-^s-isr-i, we find 

ecR{E) ^ ecR{>E) ~ (1.3 10~^eVcm-3)x£;E^V' 
where ii^EcV is in units of 10^^ eV. Thence: 



Boq — 2.2nGi?gjv^. 



(4) 



It should be stressed that such a magnetic 
field is much larger than originally assumed for 
UHECR propagation, and also larger than what 
is obtained with current models of extragalactic 
magnetic field generation (e.g. |SE|)- Therefore, 
if wave amplification by cosmic rays in the ex- 
tragalactic medium is indeed an efficient process, 
it may be the main source of magnetic field in 
the universe (outside galaxies and clusters), which 
would be very interesting in itself and would also 
have important consequences on the global field 
structure (wave content and topology). 

Clearly, the detailed study of this process 
should involve an analysis of the wave growth 
away from the linear regime as well as the non- 
linear decay (through a turbulence-like cascade 
in wavenumbers) and the isotropization of the 
wave system (which is required since the parti- 
cles are defiected by waves perpendicular to their 
stream, while they grow waves parallel to it). 
Such a study is postponed to a forthcoming work. 
We note, however, that in the above simplified 
approach, the growth timescale of the magnetic 
field at a scale resonant with extragalactic CRs 
is smaller than the age of the universe. Normal- 
izing the extragalactic density, tieg, to the bary- 
onic density of the universe, Tib, such that Jlb/i^ = 
0.02, we have wa — 4 10^ -BnG("-EG/"-b)~^^^ cm/s. 
Besides, 
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For a resonant energy -Ercsi Eq. then gives 
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Note that tieg may actually be much smaller than 
lO^^nb, because of the huge overdensity in the 
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Figure 1. Angular spread of protons of 10^® eV (left) and 10^° eV (right), propagating in a 10 nG 
turbulent magnetic field with Ac = 1 Mpc, after different times of flight, as indicated. 



galaxy clusters, so that magnetic fields could in- 
deed have been built, by now, by the CR popula- 
tion itself, up to scales resonant with very high en- 
ergies. For completeness (although we stress once 
more that the above is a very crude approach), we 
give the resonant scale for a CR of energy E in 
the self-consistent field Bcq{E) (Eq. 0}: 

Are. ^0.5 Mpc ^ . (7) 

3. Particle transport in turbulent fields 

Before turning to the possible influence of mag- 
netic fields on the UHECR spectrum, we consider 
the general problem of CR transport in turbulent 
magnetic fields. In a galaxy like ours, the turbu- 
lent component is roughly as intense as the regu- 
lar one (or even higher), with values around 5/iG, 
and a wavelength distribution compatible with a 
Kolmogorov spectrum, up to scales of 50-100 pc. 
In galaxy clusters, the observational data is very 
scarce, but SB /B may also be of order 1, with co- 
herence lengths of 10-100 kpc, over scales of a few 
Mpc. As for the deep extratragalactic space, cur- 
rent limits based on Faraday rotation measure- 
ments depend on the coherence length of the field, 
with possible values spanning a very wide range 
from totally negligible to a few tens of nC |7I8| . 
In all the models considered, the regular magnetic 
field in extragalactic voids is expected to be very 
low, and thus if magnetic fields are high enough 



to influence the UHECR spectrum, they should 
be dominated by the turbulent component. 

The transport of charged particles in a tur- 
bulent field is conceptually very simple, since it 
comes down to integrating trajectories influenced 
by the sole Lorentz force. However, its detailed 
treatment is limited by the nature of the prob- 
lem, which is essentially chaotic (and of course 
by our ignorance of the real configuration of the 
fields) . We are thus forced to keep to a statistical 
description of the trajectories. 

One should first note that individual trajecto- 
ries only depend on the particle gyroradius (or 
energy per charge). In weak fields or at high en- 
ergy (cf. Eq.El, when the particles have gyroradii 
much larger than the coherence length, Ac, they 
are only slightly deflected over that length, with 
S9{Xc) ~ Tg/Ac- Since the direction of the deflec- 
tions are not correlated from one coherence length 
to the next, particles diffuse in angle (relative 
to the initial direction), with an angular spread 
ag ^ t^/^. Figure n shows that Gaussian profiles 
fit very well the proton angular distributions at 
various times, as obtained with a numerical sim- 
ulation. It also appears that the angular decor- 
relation occurs on a longer timescale for higher 
energy particles, because of their larger rigidity. 
Figure El also shows that the isotropization pro- 
cess spreads over typically one order of magnitude 
in energy, and two orders of magnitude in time. 

Once the particles are completely isotropized, 
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Figure 2. Angular correlation, {cos 59), of a set of protons propagating in a 10 nG turbulent field, as a 
function of energy, at different times (left) and as a function of time, at different energies (right). 



it is appropriate to abandon the description in 
terms of individual trajectories, and consider 
their propagation as a diffusion process, in which 
the average linear distance traveled by CRs in- 
creases as Ar^ — 6DAt (this equation actually 
defines the diffusion coefficient, D). To analyze 
in greater detail the transition from a ballistic 
to a diffusive propagation regime, we have per- 
formed a numerical simulation of the trajectory 
of charged particles in a purely turbulent field, 
represented by a sum of Nm modes as [0]: 

B = ^ Ak^in exp{iknz'^ + if3„), (8) 
1 

where f„ = cosa„xJj + «sina„y^, a„ and /3„ 
are random phases (chosen once for all) and 
[x'^, y'„, z'J = [TZ{On,(f>n)] X [x, y, z] are coordinates 
obtained by a rotation of the reference frame 
bringing the z axis in the direction of the n*'' 
contributing wave (i.e. k„ is in direction [On, 
also chosen randomly). The amplitude A{kn) is 
determined as a function of ||k„|| according to 
a specific assumption on the type of turbulence 
(here, we assume a Kolmogorov spectrum with 
wavelengths between Amin and Amax)- 

Then, in exactly the same spirit as in ref. 110], 
we calculated the average quantity Ar^/6At - 
which may be called the instantaneous effective 



diffusion coefficient (lEDC) - for 5000 trajecto- 
ries, and plot the result as a function of elapsed 
time, At. Results are shown on Fig.Ofor various 
levels of precision in the field description (num- 
bers of modes per decade of wavenumbers, and 
ranges of wavelengths considered). The lEDC 
first increases linearly with time (as CRs prop- 
agate in straight line with velocity c) , and then a 
transition occurs to a diffusive regime where the 
lEDC is constant, with a value which can be iden- 
tified with the diffusion coefficient, D{E), of the 
order of 2 Mpc/Myr in that case. 

It is important to note that CR propagation 
cannot be described accurately without a suffi- 
ciently precise definition of the magnetic field. 
With only 10 modes per decade of wavenum- 
bers, the transition towards diffusion occurs much 
later, and the obtained diffusion coefficient is 
much larger (both by an order of magnitude, see 
Fig. El, than if one uses 100 modes per decade 
or more. Our calculations suggest that ~ 100 
modes/dec offers a good compromise between ac- 
curacy and computation time. This casts some 
doubt on results obtained with propagation codes 
using a field described as constant over cubic cells 
of size Ac, as often found in the literature. 

Of course, the time at which the diffusion 
regime settles and the value of D depend on the 
CR energy: lower energy particles enter earlier 
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Figure 3. Instantaneous effective diffusion coeffi- 
cient (lEDC), Ar^/6Ai, as a function of elapsed 
time, for a 10^^ eV proton in a 10 nG field de- 
scribed with various numbers of modes per decade 
of wavelengths (in the indicated range). 



in the diffusion regime and have a smaller diffu- 
sion coefficient (see Fig. 0] and next section). The 
diffusion regime can usually be modeled as a ran- 
dom walk at constant velocity (here v = c) with 
an isotropic redistribution of the direction after a 
"mean free time" , which can be identified to the 
scattering time, Tg, after which the direction of 
the particle is (on average) dccorrclatcd from the 
initial direction (cf. Fig. I^J. The diffusion coeffi- 
cient can then be written as £> ~ ^'-^'''s — ^^sC, so 
that the energy dependences of D and Tg are the 
same. Note that the standard Bohm diffusion co- 
efficient is obtained for the smallest possible mean 
free path, As = rgi L'Bohm(-E) = icrg(i?). 



4. Diffusion coefficients 

In Fig. 01 we show the evolution of lEDCs for 
different particle energies. A change of behaviour 
can be observed when the energy decreases. For 
the field parameters used here (Ac — 1 Mpc, 
most probable value of B = 10 nG, or {B) ~ 
11.3 nG), a bump appears in the curves for par- 
ticles of energy lower than £'o ~ 3 10^^ eV, cor- 
responding to Tg ~ Ac/27r (Ac = Amax/5 for a 
Kolmogorov spectrum ll ). For energies larger 
than Eo , the diffusion process settles as described 
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Figure 4. lEDC for particles of various energies, 
from IQie to 10^° eV, in a 10 nG field with a 
maximum wavelength Amax = 1 Mpc. 



above, from coherent deflection to angular dif- 
fusion, isotropization and spatial diffusion. At 
lower energy, the gyroradius is small and the dif- 
fusion process can settle within a distance smaller 
than Ac. As resonant interactions occur with 
magnetic modes having wavelengths A ^ Tg < Ac, 
the particles start to diffuse in a region (or on a 
scale) where the average field does not vanish. A 
one-dimensional diffusion then settles along the 
residual mean field, with _D|| > D±, and the 
full 3D diffusion is reached on larger timescales, 
when particles explore different uncorrelated cells 
of magnetic field. This causes the bumps seen on 
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Fig.Eb- 

Quantitatively, the diffusion coefficients are 
given by the values of the plateau in Fig. ^ 
They are gathered in Fig. showing the smooth 
transition between the quasilinear regime, where 
D{E) cx £'1/3, and the non-resonant regime, 
where D{E) oc E'^, in excellent agreement with 
theoretical expectations (e.g. jjl]). The low- 
energy limit is well described by the approached 
formula D{E) ^ L'Bohm(-Bo) x {E/Eq)^/^, where 
rg{Eo) = Ac/27r, and the high-energy limit by 
D{E) = i?Bohm(£i) X {E/EiYl\ where rg(£i) = 
Ac/4. The sum of these two functions provides a 
reasonable fit of D{E) at all energies (see Fig.|S^). 

An even better fit is obtained by adding a 
Bohm-like component, D{E) — D^{E/E^), and 
with only one critical energy, 

1/3 / TTi \ / \ 2' 



D{E) = 



where 

r^(E,) 



Ac 
5 



and £>, = -crg(£',). 



(9) 



(10) 



The accuracy is better than ~ 10% (see Fig.^t). 

As can be seen, the diffusion coefficients can 
be obtained straightforwardly from one single pa- 
rameter, which gathers the relevant informa- 
tion about the magnetic field intensity and coher- 
ence length. 




E (in eV) 



Figure 6. Decorrelation time and scattering time 
as a function of energy: io.9j ii/2 a-nd io.i ai'e the 
time such that (cos(6'(t) - 6'(0))) = 0.9, 0.5 and 
0.1, resp.; = iD/c^ (see text). 



As the slope of D{E) in logarithmic scales goes 
from 1/3 to 2, one may identify an approximate 
Bohm scaling (slope 1) for a limited range of en- 
ergies TU!. However, even there the Bohm diffu- 



sion coefficient, gTgC, is never reached, the "clos- 
est approach" being for Eq, with D ~ SDeohm- 
It should be noted here that while Casse et al.|10| 
found D{E) cx E"^/^ for E Eq, our more 
precise calculations confirm the theoretical ex- 



7 




1 10 100 1000 10' 10^ 10= 10' 



At (Myr) At (Myr) 

Figure 7. Left: instantaneous effective diffusion coefficient (same as Fig. taking into account proton 
energy losses. Right: Distance traveled by CRs of different energies, with and without energy losees, as 
a function of time. 



pectations in the non-resonant regime, namely 
D{E) cx . This scaling is also confirmed to 
be due to the E"^ dependence of the scattering 
time. In Fig. (HJ we plot the angular decorrela- 
tion times, ti/2 and io.i, derived from Fig. [3 
such that {cos 59{t)) = 0.9, 0.5 and 0.1, respec- 
tively. Also plotted is the scattering time, defined 
from the random walk modeling of the diffusion 
process: = ?>D/c^. As can be seen, Ti/2 gives 
an excellent definition of Tg. An E^ law is also 
shown to perfectly fit the curves. 

5. Transition from ballistic to diffusive 
regime 

As shown above, the propagation of CRs in a 
turbulent magnetic field smoothly passes from a 
ballistic regime, where the deflections are small 
and the distance traveled by the CRs away from 
their sources grows as r ~ ct, to a diffusive regime 
where it goes (on average) as r ~ V^Di. In prin- 
ciple, the diffusion regime can always be reached, 
provided that one waits for a long enough time 
(a few D/c^). To be specific, one may define 
Tdiff(-E') = AD{E)/c^ as the time required to 
reach the diffusion regime. Likewise, Adiff(-E) = 
CTdiff = V4£'Tdiff' — AD{E)/c is the distance trav- 
eled away from the source before CRs of energy 
E isotropize and diffuse. 



However, CRs may not "survive" long enough 
to enter the diffusion regime. When energy losses 
are included in the propagation code, the lEDC 
curves look like in Fig. |7^. A diffusion plateau 
is never reached, because the particles lose en- 
ergy and their instantaneous diffusion coefficient 
then drops to lower and lower values. In Fig. da, 
the traveled distance is shown to "saturate" , as 
CRs travel a negligible distance at low energy. 
With the assumed field parameters, all particles 
above 10^" eV have the same transport proper- 
ties, because they pass below the photo-pion pro- 
duction threshold while they are still in the bal- 
listic regime. 

Interestingly, for the field parameters consid- 
ered here, the transition between ballistic and dif- 
fusive regimes for particles in the GZK range (i.e. 
with energies between 3 10^^ eV and 10^*^ eV), oc- 
curs just around the GZK horizon, Rqzk{E) — 
CT\oss{E) (where the energy loss time scale is de- 
fined as Tioss(^') = E/E). This is illustrated in 
Fig-ISt: where we plot Rgzk{E) for protons, ''He, 
-'^O and ^^Fe nuclei, together with the above- 
defined diffusion length, Xdis{E). This change 
of propagation regime in the GZK region of the 
CR spectrum is one of the ingredients of a possi- 
ble alteration of the standard GZK feature, as we 
now discuss (see for more details). 
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6. Magnetic horizons and spectrum modi- 
fication 

To see how magnetic fields can modify tfie ap- 
pearance of the CR spectrum around the GZK 
feature, usually calculated with B = 0, it is use- 
ful to introduce the concept of magnetic horizon, 
i.e. the maximum distance which (most) CRs can 
travel away from their source in a given magnetic 
field. This is set by the energy loss time, rioss 
(or the age of the source, ts, if it is smaller), and 
the diffusion coefficient (or the lEDC if the diffu- 
sion regime is not reached). For high energy par- 
ticles, roughly propagating in straight line, the 
magnetic field has no influence and one simply 
gets the usual GZK horizon. For low energy par- 
ticles, the diffusion process prevents significant 
propagation beyond the magnetic horizon: 

(E) ~ ^4D{E)noss{E). (11) 

According to the standard GZK argument, the 
CR flux above 10^° eV should be sharply re- 
duced compared to that below 3 10^^ eV, say, be- 
cause of the sudden decrease in the GZK horizon 
(cf. Fig. |SJi). However, in the presence of rela- 
tively high magnetic fields, low-energy particles 
also have a limited range. In this sense, mag- 
netic fields behave as "low-cut filters", whereas 
the CMB makes the universe a "high-cut filter" 
(GZK effect). The association and tuning of such 



filters (in series or in parallel) can lead to various 
shapes of the propagated CR spectrum, different 
from the simple, "universal" GZK feature. 

In Fig.lHb, we plot both the magnetic and GZK 
horizons for protons, as a function of energy. Of 
course, Eq. Hll() does not apply when rioss < ''diff, 
equivalent to -Rmagn > Rgzk, so that the actual 
CR horizon is always the smallest of the two val- 
ues (lowest curve on Fig. IHb). In the diffusion 
regime, one may rewrite the horizon radius as 
m = Adiff(Tioss/Tdiff)^/^ or numerically: 

rH 0.58 Mpc ^ ( ^] Y (12) 
ZB^r^ VlMyr; ^IMpc^ ^ ' 

As can be seen on Fig. IHt, in a 10 nC field, 
CRs of 10^^ eV cannot come from sources more 
distant than 100 Mpc, i.e. a volume 1000 times 
smaller than in the _B ~ case. The same num- 
ber of sources should thus contributes at 10^^ eV 
and 10^" eV! This does not mean, however, that 
there should not be any GZK feature in that 
case. Indeed, while magnetic fields prevent CRs 
from diffusing far from their sources, they also 
increase the CR density around each source, i.e. 
within the magnetic horizon. This is nothing but 
the magnetic confinement effect, very familiar for 
Galactic CRs. Since the number of particles re- 
mains the same, both effects exactly compensate, 
provided that the magnetic confinement spheres 
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Figure 9. Filtering of low-energy CRs by magne- 
tized filaments (left) and out of the SGP (right) . 



(of radius th) centered on the different sources 
are big enough to merge. If this is the case, the 
propagated CR spectrum with and without mag- 
netic field are identical: the spatial horizon of the 
original GZK argument simply translates into a 
time horizon, CRs at 10^^ eV being able to con- 
tribute for a much longer time to the observed 
flux than CRs above 10^° eV. In fact, the propa- 
gated spectrum only depends on the time-of-flight 
distribution of the detected CRs, and it obviously 
makes no difference whether the trajectories were 
curved or not (see also ref. |T^). 

If the magnetic confinement spheres do not 
fully merge, however, interesting spectral effects 
can appear. The exact shape of the UHECR spec- 
trum thus depends crucially on the source gran- 
ularity, i.e. the typical distance between most 
important sources, ARg [Oj . 

We simply mention here a few typical effects. 
First, if we are outside the confinement sphere 
around a source, it is clear that we cannot re- 
ceive any CR from it. This may be of cru- 
cial importance to understand the transition be- 
tween Galactic and extragalactic CRs. Accord- 
ing to some models this could occur around 
3 10^^ eV. However, CRs of that energy have a rel- 
atively low diffusion coefficient, and thus a small 
magnetic horizon. In order for the source dis- 
tance, Ds, to be smaller than rn, the magnetic 
field (between the source and the detector) must 
be lower than (see Eq. I12|) : 
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(13) 



which is roughly the "equipartition field" at that 
energy (Eq.m. For such values, 310" eV CRs 
cannot come from sources more distant than 
5 Mpc (or 3 Mpc if S = 10 nC), which leaves 
very few potential sources. This condition can 
be made even stronger if one considers that the 
source and the detector are in much larger fields. 
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Figure 10. CR confinement in magnetized struc- 
tures. Top: GCR halo (5 kpc size), with different 
values of the field coherence length. Bottom: ex- 
tended halo (assumed size of 2 Mpc, e.g. a joint 
superhalo with M31), with Ac — 10 kpc and dif- 
ferent values of the superhalo magnetic field. 



namely a few /iG inside galaxies. From a gen- 
eral point of view, CRs first have to get out of 
their local environment, then diffuse towards our 
Galaxy, and finally reach the Earth - all within 
Tioss (or ts if it is smaller). Trajectories in high 
field regions are longer, which tends to suppress 
the low-energy CRs (low-cut filter). 

Other interesting situations may arise from 
sources located behind a filament or a galaxy clus- 
ter with high magnetic field. As illustrated on 
Fig. high energy particles are not affected, 
while low-energy CRs get a much longer path (or 
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can even be reflected), which increases their en- 
ergy losses. This will smoothen the GZK suppres- 
sion, and in some cases prevent the lowest energy 
particles to reach us at all. By observing the CR 
spectrum coming from behind a highly magne- 
tized filament or cluster with sufficient statistics, 
one should be able to observe a magnetic shadow 
(in the form of a modified spectrum). 

An other interesting effect appears if most 
sources and our Galaxy are inside a contiguous 
region of relatively high magnetic field (e.g. the 
supergalactic plane, SGP ^Sj)- Low-energy CRs 
then diffuse out of the SGP on a scale comparable 
to its thickness, while high-energy CRs can reach 
us (cf. Fig. |5Jd). Magnetic fields thus again be- 
have as low-cut filters, modifying the propagated 
spectrum. 

Finally, we note that the magnetic confinement 
of GCRs generally applies to particles with a 
diffusion length smaller than the height scale of 
the so-called cosmic-ray halo, inferred from low- 
energy CR phenomenology to be of the order of 
~ 5 kpcfn]. It is shown in Fig. llOh . that this cor- 
responds to CRs of energy lower than ^ 3 10^^ eV, 
i.e. up to the ankle of the CR spectrum. How- 
ever, one may consider an extended halo with a 
much smaller magnetic field but a larger scale, 
where CRs of higher energy could still be con- 
fined, although with a smaller enhancement fac- 
tor. In Fig. llOb . we show that a magnetic field 
of 300 nG could confine protons in an extended 
halo of ~ 2 Mpc up to GZK energies, which means 
that Galactic sources might be relevant even for 
the highest part of the spectrum (of course, col- 
limated or intermittent sources, such as GRBs, 
would then be required, in order to avoid a di- 
rect irradiation of the Earth, which would violate 
current limits on the CR anisotropy at these en- 
ergies). This is all the more true if the most en- 
ergetic CRs are heavy nuclei. An extended halo 
with a field of a few tens of nG is then sufficient 
to confine UHECRs injected by Galactic sources. 

In conclusion, the study of the UHECR spec- 
trum should involve a discussion of the magnetic 
field inside and outside galaxies. 
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